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are large regions consistent with current observable bounds. Whenever the interac- 
tions dominate, we discover significant deviations from the relations /nl ~ ''"dec 
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nated by (?NL rather than by /nl- Due to self-interactions, both the relative signs of 
/nl and ^^nl and the functional relation between them is typically modified from the 
quadratic case, offering a possible experimental test of the curvaton interactions. 
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1. Introduction 

In the curvaton scenario [Q, primordial perturbations originate from quantum fluc- 
tuations of a curvaton fleld a which remains effectively massless during inflation 
and contributes very little to the total energy density. After the end of inflation, 
the curvaton energy density stays nearly constant until the fleld becomes effectively 
massive, while the dominating radiation component generated at reheating scales as 
Prad oc a~^. The curvaton contribution to the total energy density therefore grows 
rapidly after the end of inflation and its perturbations start to affect the expan- 
sion history. In this way the initial isocurvature fluctuations of the curvaton fleld 
get gradually converted into adiabatic curvature perturbations. The observed pri- 
mordial perturbation can originate solely from the curvaton perturbations, although 
scenarios with mixed inflaton and curvaton perturbations are also possible 0. As 
the curvaton flnally decays and the decay products thermalize, the adiabatic hot big 
bang epoch is recovered and the curvature perturbation freezes to a constant value 
on superhorizon scales. 

Predictions of the curvaton scenario depend crucially on the form of the curvaton 
potential 0, |, |, g % |, |, H (and the background evolution |2[). Although 
the curvaton must be weakly coupled to the thermal bath after the end of inflation, 
interactions of some type should be present in realistic models. 
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In the present paper, we consider self-interacting curvaton models with the 
generic potential 

V = Im'a^ + Aa"+^ , (1.1) 

with n = 0,2,4,6 and A > 0. We set Mp = (SttG)"^/^ = 1 throughout the paper. 
In [|l^ it was found numerically that the predictions of such models can deviate 
significantly from the extensively studied quadratic case [|^, as shown already in 
[^, in the limit of small interactions. For non-quadratic potentials the amplification 
of curvaton energy after inflation is less efficient than in the quadratic case, but 
one may still generate the observed amplitude of primordial perturbation since the 
curvaton perturbation generated during inflation, Sa^/a^, can be larger than 10~^. 
As shown in [|10[, the correct amplitude can be achieved in a relatively large part of 
the parameter space. 

However, the dynamics is very complicated and, for non-renormalizable poten- 
tials, the curvature perturbation ( = AN (here understood to contain all orders of 
perturbation theory) displays oscillations as a function of the initial curvaton field 
value (T*. This refiects the dynamics of transition from curvaton oscillations in the 
non-quadratic part of the potential to the quadratic part, as discussed in |10[. For 
the (marginally) renormalizable case, n = 0, the transition does not give rise to an 
oscillatory behaviour of ( but nevertheless affects the final value of the perturbation. 

Here we extend the analysis of [0 and focus on primordial non-Gaussianities 
generated in the class of curvaton models with the potential ( |1.1|) . We use the AN 
formalism and solve the equations of motion numerically. We perform a system- 
atic scan through the parameter space and evaluate the non-linearity parameters 
/nl and (7nl- Combining this information with the amplitude of perturbations, we 
find the regions of the parameter space that are consistent with current observa- 
tional constraints. As a result of the oscillatory behaviour of C(c"*)) the observational 
constraints on the self-interacting model ( [L.ID differ significantly from the quadratic 
case. 

The paper is organized as follows. In Sect. ^ we give our basic definitions and 
discuss the generic features of the non-linearity parameters /nl and ^^nl- In Sect. |], 
we tackle the special case of the renormalizable, four point interaction, for which we 



can use the analytical estimates derived in [^. In Sect. Q we extend our analysis 
to non-renormalizable interactions with n = 2, 4, 6 by resorting to numerics. This 
section contains our main results. Finally, in Sect. B we present our conclusions. 



2. Non-Gaussianities in self-interacting curvaton models 



We use the 6N formalism [0, |T5| and assume that the curvature perturbation arises 
solely from perturbations of a single curvaton field generated during infiation. The 
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curvature perturbation ^ can then be expanded as 

at, x) = N'{t, t,)5a,ix) + lN"{t, t,)5a,{xf + \N"'{t, t,)5a,{xf ■■■ . (2.1) 

2 D 

Here N{t,t^) is the number of e-foldings from an initial spatially flat hypersurface 
with fixed scale factor a{t^) to a final hypersurface with fixed energy density p(t), 
evaluated using the FRW background equations. The final time t should be chosen as 
some (arbitrary) time event after the curvaton decay. The prime denotes a derivative 
with respect to the initial curvaton value cr*. Here we take t^ to be some time during 
inflation soon after all the cosmologically relevant modes have exited the horizon and 
assume the curvaton perturbations 6a^ are Gaussian at this point. The expansion 
( p.l|) is then of the form 

at, x) = Cg(t, x) + ^hLQ{t, xf + ^g^LQit, xf + ---. (2.2) 

where Cg{t, x) is a Gaussian field and the non-linearity parameters are given by 

5 A^" 

•/'nl = 6 iV^ (2.3) 
^NL = . (2.4) 



Here we neglect all the scale dependence of the non-linearity parameters |]I6|. With 
this assumption and neglecting higher order perturbative corrections, the constants 
/nl and (7nl measure the amplitudes of the three- and four-point correlators of ( 
respectively. 

We assume the curvaton obeys slow roll djTiamics during inflation and introduce 
a parameter r* to measure its contribution to the total energy density p at 



Pa 

r* = — 
P 



The inflationary scale if* appears as a free parameter in our analysis, up to certain 
model dependent consistency conditions. Assuming inflation is driven by a slowly 
rolling inflaton field, we need to require ^ lO^^^e in order to make the inflaton 
contribution to ( negligible. In this setup we also need to adjust e = —H^jHl, 
determined by the inflaton dynamics, to give the correct spectral index, n — 1 = 
2e — 2r]aa [|^. The curvaton contribution, r/o-o- = V"{<^*)/3H^, is typically negligible 
because of the subdominance of the curvaton. The curvaton mass needs to be small 
but the same also holds for the inflaton mass. We assume this minimal setup in 
the current work for definiteness since our main goal is to discuss the new features 
arising from curvaton self-interactions. 

After the end of inflation, we assume the inflaton decays completely into radiation 
and the universe becomes radiation dominated. We introduce a phenomenological 
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decay constant T to account for the coupling between the radiation and the curvaton 
component. The evolution of the coupled system is given by 



a + (3if + r)(T + mV + A(n + 4)(t"+^ = (2.6) 
= -4:Hp, + Tct^ (2.7) 

3H^ =Pr + Pa . (2.8) 

The initial conditions are given by pr = 3H^ and = l^(cr*) = r*/(l — r*)/?!- 
specified at time corresponding to the end of inflation. We also put & = 0. Given 
the parameters n, A and m, which determine the potential ( p.. ID , and the two initial 
conditions r^, and H^,, we can calculate in ( ^.11 ) from this set of equations. To 
find the curvature perturbation, we set 6a^ = H^/{2it) and compute ( = N{a^ + 
(5cr*) — A^(cr*). For a given set of parameters, we adjust the decay width T so that 
the observed amplitude is obtained, ( ~ 10^' 



We treat F as a free parameter since we have not specified the curvaton couplings 
to other matter, in particular to the Standard Model fields. However, since the 
primordial perturbations have been observed to be adiabatic to a high degree [0, 
the curvaton should decay before dark matter decouples in order not to produce 
isocurvature modes. Assuming the freeze-out temperature for an LSP type dark 
matter model with Tdm ~ (9(10) GeV, this translates to a rough bound 

F > lO^^^GeV = 10"^^ . (2.9) 

While this bound could be relaxed in non-minimal constructions, we assume it for 
definiteness for the rest of our work. 

2.1 Analytical considerations 



In Sect. § we solve the set of equations (|2.6| ) - ( p.8|) numerically and compute the 



resulting curvature perturbation (p.l|) and the non-linearity parameters (|2.3|) and 
( p.4|) . However to gain some physical insight of the results, it is useful to start by 
discussing generic approximative analytic expressions that can be derived for /nl and 
S^NL- Assuming the curvaton decays instantaneously at H^ec = F and neglecting 



the coupling between curvaton and radiation before t^ec, one can can estimate the 
non-linearity parameters by H, 



/nl = ( 1 + ^f^) - I - ^ (2.10) 



_ 25 



^i^ + %£^)-ilfl + ^^)+ (2.11) 



^' '^;JV 4 + 16 
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Here rdec = pa/ p\dcc ~ cr^gj.(m/r)^/^ and Uosc is the envelope of the oscillating curva- 
ton field at some time tosc after the beginning of oscillations in the quadratic part of 
the potential. The results are independent of the precise choice of tosc- 

If the curvaton potential is exactly quadratic, cTosc cr* and the non-linearity 
parameters /nl and ^^nl are uniquely determined by rdec, i-e. by the curvaton energy 
density at the time of decay. However, interactions in general make the function 
o"osc(o"*) non-hnear and, especially in the limit r^^c ^ 1, the quadratic predictions can 
be greatly altered due to the derivative terms in (|2.10| ) and (|2.11|) . If the interactions 



dominate the curvaton dynamics at the time of inflation, the non-linearity parameters 
can only very crudely be approximated by /nl ~ ''"dec ^^^1 (7nl ~ "'"dec- ^^^^ t)e 
discussed in Sect. ^, in the interacting case the non-linearity parameters can deviate 
significantly from these naive estimates. This follows from the non-trivial behaviour 
of derivatives of o"osc(o"*) in (|2.1CI|) and 



In particular, we find that for n = 2, 4 in the potential ( |1.1| ) the signs of the 
non-linearity parameters /nl and (^nl oscillate as a function of a*. Therefore, even if 
'"dec <^ 1, it is always possible to find regions where the non-Gaussianities are acci- 
dentally suppressed and do not confiict current observational bounds. For example, 
from ( |2.10|) we see that /nl ~ for rdec ^ 1 whenever foscC"osc/(^osc)^ ~ — 1. Using 
the definitions (p. 31) and (|2^) we find, 



/nl(«^*) = V (^(7NL - ^/^l) , (2.12) 

which gives an estimate |A(t| ~ {g^i,N')~^ for the typical size of the regions /nl ~ 0. 
This is larger than the range of field values probed by the curvaton fiuctuations 
produced during infiation, \5a^:\ ~ if g^i, < 10^. Therefore, for g^i^ < 10^ 

we may conclude that the regions /nl ~ correspond to a self-consistent choice 
of initial conditions and are not destabihzed by the fiuctuations produced during 
infiation. 



3. Renormalizable potential 

The special case of a quartic interaction term in the potential, 

1 



V = ^mV^ + Xa^ , (3.1) 



can be discussed using the analytical estimates presented in [|I0|. If the interaction 
dominates over the quadratic term at the time of infiation, the curvaton oscillations 
start in an effectively quartic potential. As the amplitude of the oscillating field 
decreases, the quartic term dilutes away and the potential becomes quadratic. As- 
suming the universe remains radiation dominated until the curvaton decay, rdec ^ 1, 
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the envelope of the oscillating curvaton in the quadratic regime is given by 



(3.2) 



The function crosc(cr*) can be estimated using Eq. (4.9) in ||T0[, which assumes Xa^ > 
and is derived to leading order precision in r^ec ^ 1. In approaching the quadratic 
limit, Acr^ ^ m^, the analytical estimates of [10] can no longer be used to obtain 
quantitative results but the qualitative features remain correct. 
To leading order in r^ec, Eqs. ( |2.1U| ) and ( [^.lip read 



NL 



5'NL 



3 '"dec 

50 

^ " dec 



// 

O^oscC'osc 



a, 



osc/ 
/// 



'^osc'^osc I 3croscO"osc 



[a, 



[a, 



(3.3) 
(3.4) 



and by substituting ( p.2|) into these, we find (semi) analytical predictions for /nl and 
(7nl- The results are illustrated in Fig. (p. 
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Figure 1: Contour plots of /nl (left panel) and ^nl (right panel) with variables r^ec and 
^/Xa^/m as x and y axes, respectively. On the left panel, the contours run from (black) 
to 100 (white) with spacing of 10. On the right panel, they run with spacing of 500 from 
-5000 (black) to (white). 

Although the n = case does not result into oscillatory behaviour of C('^*); 
non-monotonous features appear when considering derivatives of ({c^:), i.e. the co- 
efficients N",N"' etc. in ( ^.ip . These features become the more pronounced the 
higher the order of the derivative A^^"*) is. This can also be observed in Fig. (|l|) 
where both /nl and gjssh display non-monotonous behaviour as a function of a* in 
the region Xa^ ~ m^. For /nl the non- monotonous features are very mild but (^nl 
shows a significantly more sensitive dependence on cx^,. 

As discussed in [|10|, for Xa^ ~ the transition from the quartic to the quadratic 
part of the potential takes place soon after the onset of oscillations when a and & 
still act as independent degrees of freedom. These quantities in general do not have 
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a similar dependence on the initial field value cr* . Variations of cr* therefore affect the 
effective equation of state of the oscillating curvaton in a non-trivial fashion and this 
is the origin of the structure seen in Fig. (|I]). If Xa^ 3> or Xa^ -C m^, no similar 
structure is seen. In the former case, the transition into quadratic region happens 
when the dynamics of the oscillating curvaton is already well described by a single 
dynamical degree of freedom, the amphtude a. In the latter case the potential is 
almost Gaussian from the beginning. 



4. Non-renormalizable potentials 

The discussion in the preceding section, together with the oscillatory behaviour of 



C presented in leads to an expectation that for non-renormalizable curvaton 
potentials the naii've estimates /nl ~ l/'^dec and (7nl ~ /nl ^e violated for a 
range in the parameter space. It is however not obvious at all how large this range 
might be. Since the analtyical approximations used in the previous section cannot 
be applied to non-renormalizable potentials, we use numerical methods to track the 
dynamics at hand. 

To calculate the values of /nl and (7nl for the non-quartic cases, we use code 



developed for |TD[ to compute the values numerically using Eqs. |2.3| and |2.4|. To 



obtain the derivatives, we calculate N for five different initial conditions separated 
by a spacing a, and then use five-point stencil to calculate the first, the second 
and the third derivative of N. We adjust the stepping a so that numerical noise is 
minimized. Furthermore, we adjust a for each pixel independently. 

4.1 Qualitative features and differences from naive expectations 

For all choices of m and n, the parameter space is divided into two areas by a line 
where the quadratic and the non-quadratic terms are equal for the initial curvaton 
field value, i.e., 

2 2 \ ^* 

-m erf = A . 

2 * M" 

Below this quadratic line we should recover the behaviour predicted analytically 
for the quadratic case. Above this line the self-interactions modify the behaviour 
substantially. 

For n = 6 case the behaviour in the non-quadratic regime is smooth and qualita- 
tively in close resemblance with the case n = 0, which was described in the previous 
section. This is due to the fact that if n > 6, the field does not oscillate in the 
non-renormalizable part of the potential at all. However, for the cases n = 2 and 
n = 4 large oscillations of C as a function of its parameters and initial conditions 



are present in the non-quadratic regime as described in |]I0[. Here the expressions 
/nl ~ l/'^dec and /^l ~ fi'NL provide only very rough estimates of the non-linearity 
parameters, and in the interaction dominated region the actual values of /nl and (^nl 
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can deviate significantly from tliese estimates. Tliis beliaviour can be traced back to 
tfie terms including derivatives of cTosc in the Eqs. (|2.1CI| ) and 
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(a) |/nl| plotted against the intial conditions (b) |/nl| plotted against rdcc for a fixed value 
and iJ,. = 5 X 10^^, corresponding to the vertical 

stripe in Fig. 2(a)| . Red points correspond to 
/nl > and blue points to /nl < 0. 
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(c) |/nl| plotted against rdcc for aU values of i?* (d) |5nl| plotted against |/nl| for all the points 
in 



2(a). Red points correspond to /nl > and in Fig. 2(a). The green line corresponds to the 

linear relation (Jnl ^ /nl and the blue line to 



blue points to /nl < 0. 

the quadratic relation ^nl ^ /nl- 

Figure 2: The behaviour of |/nl| and I^'nlI for n = 4 and m = 10^^^ 



In Fig. |2(a)| we have plotted the value of |/nl| against the initial conditions 
and H^,. The vertical line in this figure corresponds to a range of parameters, corre- 
sponding to the fixed value of if,,, = 5 x 10~^, for which |/nl| is plotted against r^ec in 
Fig. |2(b)| . Here /nl can be seen oscillating around the 1/rdec-estimate. Oscillations 
arise from the derivative terms present in Eq. (|2.10|) , which describe the impact of 
the self-interaction terms on the dynamics of the curvaton. Most points give rise 
to a larger /nl than one would expect from the estimate 1/rdec, however since /nl 
actually changes sign, points can always be found where /nl is arbitrarily close to 
zero. In Fig. P(b)| several different values of |/nl| correspond to a given point r^ec 
since different choices of initial conditions can be degenerate yielding the same r^ec- 
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In Fig. P(c)| we plot again |/nl| against r^ec, but we no longer constrain if*, but 



instead plot this for all points in Fig. p(a)| . Again each point corresponds to a set of 



parameters producing the observed final amplitude of the primordial perturbations. 
As we allow H^, to take different values, a family of curves is drawn, where each curve 
is similar to the curve present in P(b)| , resulting into the noisy scatter present in the 



Fig. \i{c)[ It is also noteworthy that for given fixed value of rjec, there are multiple 
sets of parameters which all give the same final amplitude for the perturbations, but 
different value of /nl (and (^nl)- 

In Fig. p(d)| the value of I^^nlI is plotted against |/nl|- Here three contribu- 
tions can be clearly distinguished: The 1/rdec-relation arising while in the quadratic 
regime, the 1 /r ^g^-relation in the non-quadratic regime that is due to self-interactions, 
and scatter around those lines due to the oscillations caused by the self-interactions. 
This scatter can be understood by considering Fig. |2(b)| where /nl can be seen os- 



cillating around the analytical estimate. Plotting (/nl against rjec would produce a 
qualitatively similar plot as Fig. ^(b)| showing the oscillations around the l/vdec and 



l/r^g^ estimates deriving from the terms in eq. (|2.10 ). However, in general /nl and 



^^NL do not oscillate with the same phase, and hence when plotting (^nl against /nl, 
scatter is created. 



In Fig. |3(a)| - |3(d)| the magnitudes and signs of /nl and (^nl are plotted when 



scanned through different initial conditions. Here again the oscillatory features can be 
clearly distinguihed as both /nl and (7nl oscillate in the regime initially dominated by 
the non-quadratic interaction. Moreover it is worth emphasizing that not only does 
the absolute magnitude of /nl and (^nl show oscillatory behaviour in this regime, 
but also their signs change along the oscillations as shown in Fig. p(b)| and |3(d) 



Futhermore the oscillations of /nl and (7nl have different periods and phases, i.e. the 
zeros of /nl and (7nl are not related in a simple fashion. 



In Eq. ( p. 121) we presented a relation between the derivative of /nl with respect 



to cr^ and the values of (^nl and /nl- If /nl = 0, this takes a particularly simple 
form, /^L = 9/5A^'(7nl, which can be clearly seen in Fig. P(c)| . 

Note that there are two different sources of non-Gaussianity: one is the subdom- 
inance of the curvaton at the time of decay rdec, while the other is just the non-linear 
evolution of curvaton perturbations, encoded in the function dose in Eqs. ( p^.lOD and 
( p. 11 ). Even if r^ec ~ 1, large non-Gaussianity can be generated by the evolution of 



the curvaton. For example, from Eq. ( p. 10]) we see that if a'^^^ — )■ 0, /nl can be very 
large even though r^ec ~ 1- This can be understood qualitatively by looking at the ex- 
pression for the linear part of the curvature perturbation, ( ~ rdcc{o'Lr/<^osc)So'* \Wi ■ 



If o-Qgj, —7- 0, we see that we need to increase 5(T* to keep ( ~ 10~^. Therefore, the 
higher order terms become significant in this limit generating large non-Gaussianities, 
just like in the limit r^^c 0. 
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(a) |/nl| plotted against and 



(b) Sign of /nl plotted against r, and , where 
the red color corresponds to positive and blue to 
negative /nl- 
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(c) |<?nl| plotted against and i/* 



(d) Sign of 5nl plotted against and iJ* , where 
red color corresponds to positive and blue to neg- 
ative (7nl- 



Figure 3: Magnitude of |/nl| and |5nl| for n = 4 and m = 10 
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4.2 Allowed regions of parameter space 

In Figs. ^ - 1^ we have plotted the points in the parameter space which are compatible 
with observations of the primordial perturbations. These points give rise to the 
observed amplitude of the perturbations while producing /nl and ^^nl which are 
within the current observable limits. 

The range of the initial conditions i/* and has been chosen so that all interest- 
ing features should be within the plots. is also bounded from above, if* < 10~^, 
to prevent the excessive production of primordial tensor modes and at least an order 
of magnitude smaller to keep the infiaton perturbations negligible. 

Limits for /nl are given by the WMAP 5-year data [|1^, —9 < /nl < HI- 
Although in |Q a more stringent constraint for /nl is given as —4 < /nl < 80, 
we conservatively use the limit provided in [|18|. Regarding the limits for (^nl, we 



require that (^nl is the range —3.5 x 10 < s'nl < 8.2 x 10 as given in Note 
that these limits have been derived assuming /nl ~ 0, which here is not the case 
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in general. However, the bounds for /nl seem to be much more constraining in our 
case, and relaxing the limits for (yfNL even by an order of magnitude would not enlarge 
the allowed area of the parameter space significantly. Thus we adopt this limit for 
reference purposes^. 
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Figure 4: A schematical illustration of the different cuts limiting the allowed area in the 
parameter space. 

In Fig. I we give an illustration of the different features limiting the allowed 
area of the parameter space. The observational limits for /nl and (7nl constrain the 
allowed area in the very subdominant regions of the parameter space, as depicted in 
Fig. ^ by the line a. Other constraints shown in Fig. ^ arise from the internal consis- 
tency of the specific curvaton scenario studied in the present paper. The bound b is 
obtained because otherwise the initial perturbations would be too small to produce 
the observed amplitude. The bound c reflects the requirement that the curvaton 
should be massless, or V" < H^, which is necessary for the generation of curvaton 
perturbations during inflation. Because of the subdominance of the curvaton, the 
realistic bound should arguably be a few orders of magnitude tighter. However, a 
change of an order of magnitude moves the actual cut by a very small amount in the 
log-log plots. 

Finally, the bound d guarantees the absence of the isocurvature modes in dark 
matter perturbations and corresponds to the limit on the curvaton decay width given 

"'^ Recently, the authors of [^5| have obtained limits on ^nl without assuming /nl = by using 
N-point probability distribution function. However, their constraint on g^i^ is similar to that of 
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in (pD. 

It should be noted, that hmits b, c and d in Fig. || were aheady present in , 
so that the non-Gaussianty hmit a is the only new limiting ingredient provided by 
the current work on the space of parameters. 
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(a) n = and m = 10 



(b) n = and m = 10 
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(c) n — 6 and m = 10 



(d) n — 6 and m = 10 



Figure 5: Dark areas correspond to the allowed areas in the parameter space with —9 < 
/nl < 111 and -3.5 x 10^ < ^nl < 8.2 x 10^ forn = and n = 6. 



As explained in the previous sections, in the regions where the quadratic part of 
the potential dominates already initially, the dynamics are essentially linear, and the 
simple results apply. As a consequence, cTosc dependence in Eqs. ( ^.10[ ) and ( ^.11| ) 
disappears. Therefore these regions are characterized by the smooth continuous 
allowed shown in Figs. ^ - 0, which can be found in the lower left are in the 

plots. The total area of the allowed region depends on the values of m and n; e.g. 
for n = A and m = 10~^ this quadratic area is much larger than for, say, n = 4 and 
m = 10"^^ 

As mentioned previously, n = and = 6 do not display extended oscillatory 
regions. Hence the plots in Fig. ^ remain smooth also in the regime where the 
self-interaction dominates. 

The cases n = 2 and n = 4 are however characterized by significant oscillatory 
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features, as discussed in the previous sections. Due to these oscillations, the allowed 
region consists of isolated patches in the interaction dominated regime. As can be 
seen in the figures, the size of these patches decreases as we decrease the (bare) mass 
m, since this tends to make the curvaton more subdominant at the time of decay, 
decreasing the viable area. 
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(a) n = 2 and to = 10 



(b) n = 2 and to = 10 
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(c) n = 2 and m — 10 
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Figure 6: Dark areas correspond to the allowed areas in the parameter space where 
-9 < /nl < 111 and -3.5 x 10^ < ^nl < 8.2 x 10^ for n = 2. 



5. Discussion 

It seems plausible that the curvaton, like any other scalar field, has some self- 
interactions. These self-interactions can have a profound effect on the dynami- 
cal evolution of the curvaton field and its perturbations, as was discussed in [|lOl 
(and in @, |^)), where we studied the amplitude of the curvature perturbation 
in self-interacting curvaton models defined by Eq. In the present paper we 

have focussed on the non-Gaussianities of the curvature perturbations by computing 
the non-linearity parameters /nl and ^fNL for all the model parameters for which 

C ~ 10-^ 
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(a) n ~ 4 and m = 10 



(b) n = 4 and m = 10 
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(c) n = 4 and m — 10 
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(d) n — A and m = 10 
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Figure 7: Dark areas correspond to the allowed areas in the parameter space where 
-9 < /nl < 111 and -3.5 x 10^ < c/nl < 8.2 x 10^ for n = 4. 



When the curvaton has some self-interactions, the non-Gaussian statistics of the 
curvature perturbation can be quite different from that produced in the simplest 
model with a quadratic potential. In the quadratic curvaton model, the magnitude 
of /nl in the limit r^cc <^ 1 is determined by the curvaton energy density at the time 
of its decay, /nl ~ l/^dcc- However, the prediction of /nl can significantly deviate 
from this simple estimate if the curvaton has non-renormalizable self-interactions. 
As seen from Fig. |2(c)| , for such models the values of /nl scatter around the estimate 
^/^dec and typically end up being slightly larger than in the quadratic case. Thus a 
very subdominant curvaton is not favoured because it tends to yield a value of /nl 
which is in excess of the observational bounds |/nl| ^ 100 [|18], p^ . 

However, it is interesting to note that even if r^ec ^ 1, there exists regions in the 
parameter space with |/nl| < 0{1). This is because the value of /nl oscillates and 
changes its sign, as is illustrated in Fig. P(a)| for n = 4. However, even if /nl < 
the non-linearity parameter for the trispectrum (^nl can be very large, as can be seen 
in Fig. P(d)| . In these regions the self-interacting curvaton scenario gives rise to a 
rather non-trivial non-Gaussian statistics characterized by a large trispectrum and a 
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vanishing bispectrum. Such a situation, discussed aheady in , appears to be rather 
generic in self-interacting curvaton models, and possible for a wide, albeit restricted, 
range of model parameters. 

Another interesting feature of the curvaton model with self-interactions is that 
large non-Gaussianities |/nl| ^ C(1) can be generated even if the curvaton domi- 
nates the energy density at the time of its decay, r^^c ~ 1- Recently in it was 
shown that for rjec ~ 1, the entropy production at the curvaton decay can leave an 
imprint on the spectrum of primordial graviational waves, which in principle could be 
observable. If the curvaton had no self-interactions, such a signal would imply that 
no large non-Gaussianity could be generated by the curvaton fluctuations. This is 
clearly not the case when the self-interactions are included, which serves to demon- 
strate thet the self-interactions can significantly affect the generic features of the 
curvaton scenario. 

Another interesting feature that is clearly visible in Fig. |2(d)| is the breakdown 
of the relation of (^nl = — (10/3)/nl which holds true for the quadratic potential. 
For the self-interacting curvaton a large number of points, each corresponding to an 
allowed set of parameter values, can be seen to fall into the region with Ig-^hl > |/nl|- 
If the interactions are small compared to the quadratic part, the linear relation ^f^L = 
— (10/3)/nl gets replaced by g^i^ ~ ~/nl El- However, when the self-interaction term 
dominates, both of the above relations can be violated, as is seen from Fig. |2(d) 



There is nevertheless a tendency for the points to be concentrated around the lines 
\9nl\ ~ I/nlI and \g^i^\ ~ 

We should also like to point out that in the quadratic case with |/nl| ^ 1, 
the signs of /nl and (^nl are respectively positive and negative. In contrast, in the 
cases studied here, the signs can well be the same. This feature provides an obvious 
departure from the quadratic case and could offer an experimental possibility to 
constrain the curvaton self-interaction. 

It is evident from the figures presented in Sect, ^that the self-interacting curva- 
ton model provides us a rich tapestry of features, constrained in a highly non-trivial 
way by the observational limits on non-Gaussianity. As we have argued here, in 
the presence of self-interactions the relative signs of /nl and (7nl and the functional 
relation between them is typically modified from the quadratic case. Thus the non- 
linearity parameters taken together, in possible conjunction of other cosmological 
observables such as tensor perturbations, may offer the best prospects for constrain- 
ing the physical properties of the curvaton. 
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